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Recent experiments on La2Cu04 suggest that indirect resonant inelastic X-ray scattering (RIXS) 
might provide a probe for transversal spin dynamics. We present in detail a systematic expansion 
of the relevant magnetic RIXS cross section by using the ultrashort core-hole lifetime (UCL) ap- 
proximation. We compute the scattering intensity and its momentum dependence in leading order 
of the UCL expansion. The scattering is due to two-magnon processes and is calculated within a 
linear spin-wave expansion of the Heisenberg spin model for this compound, including longer range 
and cyclic spin interactions. We observe that the latter terms in the Hamiltonian enhance the first 
moment of the spectrum if they strengthen the antiferromagnetic ordering. The theoretical spectra 
agree very well with experimental data, including the observation that scattering intensity vanishes 
for the transferred momenta q — (0,0) and q — (7r,7r). We show that at finite temperature there 
is an additional single-magnon contribution to the scattering with a spectral weight proportional 
to T^. We also compute the leading corrections to the UCL approximation and find them to be 
small, putting the UCL results on a solid basis. All this univocally points to the conclusion that 
the observed low temperature RIXS intensity in La2Cu04 is due to two-magnon scattering. 



I. INTRODUCTION 

Indirect Resonant Inelastic X-ray Scattering (RIXS) 
is rapidly establishing itself as a new probe of elec- 
tronic excitations in solids. The recent increase in bril- 
liance of synchrotron radiation has made it possible to 
obser ve second order scattering processes as indirect 
RIXg i | ^ l ^ | 4 | 5 | B | 7 | « l » | iu | ii | i^ | ia | iM Ilfel. Moreover, the im- 

provements in the instrumental resolution (100 meV is 
achieved) allow for lower energy scales to be detected, 
making this technique in principle a powerful instrument 
to probe the low-lyi ng ele mentary excitations of solids, 
for instance magnonJ^ZEH 

In indirect RIXS, the energy of the incoming photons 
is tuned to match a resonant edge of an atomic transi- 
tion in the particular system that one sets out to inves- 
tigate. This resonance corresponds to exciting a core 
electron to an outer shell. The if -edge of transition 
metal ions is particularly useful since it promotes a Is 
core electron to an outer 4p shell, which is well above the 
Fermi level, so that the X-rays do not cause direct transi- 
tions of the Is electron into the low est 3c?-like conduction 

]3j^mj j2l3l4l5l6l7l8l9ll0lllll2ll3ll4ll5ll6l 

Due to the large energy involved (~5-10 keV), the 
core-hole is ultrashortlived and it induces an almost 
delta function-like potential (in time) on the valence elec- 
trons^^ ^® Consequently, elementary excitations of the 
valence electrons will screen the local potential, but have 
litUe time to do so. When the core-hole decays, the sys- 
tem can be left behind in an excited state. By observing 
the energy and momentum of the outgoing photon, one 
probes the elementary excitations of the valence electrons 



including, in particular, their momentum dependence. 

In the last few years, considerable theoret ical prog ress 
has been made to comprehend RIXS spectraP'^ESHIl and 
particularly in the understanding of the c orrelation func- 
tions that are measured by indirect JUX^^SMMMl^ jt is 
by now well established that indirect RIXS detects the 
momentum dependence of charge excitations that are re- 
lated to the electrons and holes in the d-shell in for in- 
stance the cuprates and manganites. Treating the scat- 
tering problem taking the ultrashort core-hole lifetime 
(UCL) into account has proved that the indirect RIXS 
intensity is proportional to the dielectric loss function 
and longitudinal spin excit ations of the electrons that 
couple to the core-hol^^MZI 

Recently, RIXS measurements performed by J. P. Hill 
and coworkers on the high-Tc cuprate superconductor 
La2_£!;Sr2;Cu04 revealed that RIXS is potentially able 
to detect transversal spin excitations -magnonJi^. The 
experiments show that the magnetic RIXS signal is 
strongest in the undoped cuprate La2Cu04. The mag- 
netic loss features are at energies well below the charge 
gap of this magnetic insulator, at energies where the 
charge response function S{q,uj) vanishes, as well as the 
longitudinal spin one -which is in fact a higher order 
charge response function. The proposed scattering mech- 
anism is a two-magn on sca ttering process in which two 
spin waves are createcP^'^. 

In a previous theoretical analysis we have shown that 
the magnetic correlation function that is measured by in- 
direct RIXS is a four-spin correlation one, probing two- 
magnon excitation^^. This makes indirect RIXS a tech- 
nique that is essentially complementary to magnetic neu- 
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tron scattering, which probes single magnon properties 
and two-spin correlations. In this paper, we present the 
theoretical framework of Ref. [THlin more detail and use it 
for an analysis of the experimental magnetic RIXS data 
on perovskite Cu02 layers of La2Cu04. 

We expand upon the previous considerations by pro- 
viding a detailed comparison between the theory and ex- 
periment, including also longer range magnetic exchange 
interactions in the theory -with values known from neu- 
tron scattering data. We develop the theory to account 
also for the effects of finite temperature, which give rise 
to a non-trivial single-magnon contribution to the RIXS 
signal. We also compare with the results of Nagao and 
IgarashP^, who recently computed the magnetic RIXS 
spectra based on the theoretical framework of Ref. [T51 
taking also some of the magnon-magnon interactions into 
account. 

The theory is developed on basis of the ultrashort core- 
hole lifetime (UCL) expansion. We compute leading or- 
der corrections to this expansion and show that they are 
small. This makes sure that the UCL approximation pro- 
vides a reliable route to analyze the indirect RIXS spec- 
tra. 

This paper is organized as follows: in section [ll] we 
obtain an expression for the cross section of the 2D S = 
1 /2 Heisenberg antiferromagnet in linear spinwave theory 
in terms of magnon creation and annihilation operators. 
In section IIIII we evaluate the cross section at T = 0. 
Section |IV| concerns the low temperature case. Next, 
the leading correction to the cross section in the UCL 
approximation is calculated. Section |VI| is devoted to 
the concluding remarks. 

II. CROSS SECTION FOR INDIRECT RIXS ON 
A HEISENBERG AFM 

Recently, J. P. Hill et al}'^ observed that RIXS on 
the high superconductor La2_2;Sr^Cu04 picks up 
transversal spin dynamics -magnons. In the undoped 
regime, the RIXS intensity turns out to be highest. 
The same feature was observed in the related compound 
Nd2Cu04. These cuprates consist of perovskite Cu02 
layers with a hole in the Cu 3d subshell. The low energy 
spin dynamics of these systems are properly described by 
a single band Hubbard model at half filling. The strong 
interactions between holes in the Cu Sd subshells drive 
these materials into the Mott insulating regime, where 
the low energy excitations are the ones of the S" = 1/2 
2D Heisenberg antiferromagnet: 

Ho^'^JijSi-Sj (1) 

with Jij « 146 meV for nearest neig hborsP. The su- 
perexchange integral Jij is determined from the virtual 
hopping processes concerning sites i and j: Jij — itjj/U. 
Here tij is the hopping amplitude and U is the Coulomb 
repulsion between two 3d electrons on the same site. 



In the antiferromagnetic groundstate, the Hamiltonian 
can be bosonized in linear spinwave theory (LSWT) 
where S'f i-^ 1/2 — aja^, S^' i— > and a] 
for i £ A {A being the sublattice with spin-up) and 
5| ^ b]b-j - 1/2, S+ b] and S' ^ bj for j G B (the 
spin-down sublattice). A Bogoliubov transformation in 
reciprocal space is necessary to diagonalize Hq: 

ttk = UkOk + Wkfc^k' (2) 
= "k^k + Wkolk (3) 

with 




and 

v^^sign{J^^),fi^l (5) 

where J^^ is the Fourier transform of those terms in 
Jij connecting a site in sublattice X to a, site in Y. For 
interactions up to third nearest neighbors we get 

Jk — J (cos akx + cos aky) (6) 
Jk'^ = =2 J' cos ak^ cos aky+ 

J" (cos 2akx + cos 2aky) (7) 

with a the lattice constant and J, J', J" the first through 
third nearest neighbor couplings. The final linear spin- 
wave Hamiltonian in terms of boson operators is 

Ho = const + ^ Ek (a^ttk + 4/?k) (8) 

k 

With ek = ^{J^^-J^'' + J^^r-{J^^f- 

Our aim is to understand how this picture changes 
when doing indirect RIXS. In RIXS, one uses X-rays to 
promote a Cu Is electron to a 4p state. For an ultrashort 
time, one creates a core-hole at a certain site which low- 
ers the Coulomb repulsion U on that site with an amount 
Uc- We assume that the core-hole potential is local, i.e. 
it acts only at the core-hole site. This approximation is 
reasonable as the Coulomb potential is certainly largest 
on the atom where the core-hole is located. Moreover, 
we can consider the potential generated by both the lo- 
calized core-hole and photo-excited electron at the same 
time. As this exciton is a neutral object, its monopole 
contribution to the potential vanishes for distances larger 
than the exciton radius. The multi-polar contributions 
that we are left with in this case are generally small and 
drop off quickly with distance. 

The strong core-hole potential in the intermediate 
state alters the superexchange processes between the 3d 
valence electrons. This causes RIXS to couple to multi- 
magnon excitations, as was first pointed out in Ref. ifSI 
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FIG. 1: In RIXS, a photon of momentum qin and energy 
tuned to the K-edge of a transition metal ion = ojios) 
creates a core-hole at a certain site. The superexchange in- 
teraction between this site and a neighboring other site is 
modified because the energy of the virtual intermediate states 
is changed. The same-site Coulomb repulsion U is lowered by 
Uc if the core-hole site contains no holes and is raised by Uc if 
there are two holes present. Summing the amplitudes for both 
processes, we obtain the modified superexchange interaction, 
see Eq. Q. 



The simplest microscopic mechanism for this coupling is 
obtained within the strong-coupling Hubbard model, in 
which the doubly occupied and empty virtu al states shift 
in energy in presence of the core -holePiaSl. Adding the 
amplitudes for the two possible processes shown in Fig.jl] 
lead to an exchange integral in presence of a core-hole on 
site i of 



U + Uc U -Uc 



J^J{l+r|) 



(9) 



where j is a site neighboring to i and rj = / {U^ — U^) . 
This enables us to write down the generic Hamiltonian 
for the intermediate stated 



as 



^^int — ^0 



(10) 



determined by the value of rj, because rjJ/T appears as 
a small parameter in this expansion. As for the cuprates 
J/r w 1/5 such corrections are small for the relevant 
possible values of rj. 

In the three-band Hubbard model that includes also 
the oxygen states, two important kinds of intermediate 
states appear: the poorly- and well-screened ones. Be- 
cause the Coulomb interaction of the core-hole with the 
valence electrons is large {Uc — 7.0 eV, compared to a 
charge transfer energy A = 3.0 eV^^, a copper hole can 
transfer to a neighboring oxygen to form a well-screened 
intermediate state. The low-energy sector now also en- 
compasses an oxygen hole, equally distributed over the 
ligands. We will show that, starting from a three band 
Hubbard model, Eq. ( 10 1 gives a proper description of 



both the well- and poorly-screened intermediate states, 
with 7] now a function of the parameters of the three band 
model. Before presenting these results we remark that 
scattering processes that scatter a well-screened state 
into a poorly-screened state or vice versa yield a large 
energy loss to. These are not important at low uj, where 
one will only observe scattering in the magnetic channel, 
not the charge one. 
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FIG. 2: Modification of the superexchange interaction in the 
well- and poorly-screened intermediate states. In the poorly- 
screened state, the core-hole potential Uc modifies the su- 
perexchange. For the well-screened state however, the copper 
3d hole on the core-hole site is transfered to a neighboring 
oxygen, and superexchange is only of order 0{tp^), indepen- 
dent of Uc- 



where Si creates a core-hole and sj annihilates one at site 
i. In the Hubbard framework one could identify the U 
with the Coulomb energy associated with two holes in a 
3(i -orbit al Ud = 8.8 eV, which together with Uc = 7.0 
eVlSlIlll leads to = 1.7; from U/Uc = 2/3, as suggested 
in Ref. 'M', one finds 77 = —0.8. 

The situation in the cuprates, however, is more com- 
plex and one needs to go beyond the single band Hubbard 
model to obtain a value of rj from microscopic considera- 
tions. We will do so by considering a three-band model in 
the strong coupling limit. However, it should be empha- 
sized that for the end result -the computed RIXS spec- 
trum in the UCL approach- rj just determines the over- 
all scale of the inelastic scattering intensity. As we will 
show higher order corrections in the UCL approach are 



The magnetic scattering processes for the poorly- 
screened state are very similar to the single band pic- 
ture: all copper ions have one hole and all oxygen ions 
are fiUed-shell. The superexchange processes are shown 
in Fig. ([2]). We consider the Anderson and Geertsma 
contributions to the superexchange^'* and find 



Vps 



1 



UdA\Up + 2A) 
2{2Ud + 2A + Up) \iUd - Uc){Uc ~ A)2 

1 , [l/A + l/{A~Uc)Y 



(Ud + C/c)A2 

which results in 77 = 
eV, Up = 6.0 eV, tpd 



2A -Uc + C/„ 



1, (11) 



—0.3 using the parameters Ud 
= 1.3 eV, A = 3.0 eV, and Uc 



7.0 
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eV23124|^ where tpd is the copper-oxygen hoping integral 
and Up the on-site Coulomb repulsion of two oxygen 
holes. 

The well-screened intermediate states have a similarly 
modified superexchange interaction, as shown in Fig. [2] 
Because of the large core-hole Coulomb interaction an 
electron from the neighboring oxygen atoms moves in to 
screen it, or, equivalently, the copper hole is transferred 
to the in-plane oxygen ions. Transfer out of the plane is 
not considered since the Cu M^-^-y-^ hole only couples to 
the in-plane oxygens. Because the Cu hole is transfered 
in the direction of one of its neighboring Cu ions, the con- 
tribution to the superexchange interaction for the well- 
screened state is of second order in tpd, instead of fourth 
order between two Cu sites (see Fig. [2). The rotational 
invariance around the core-hole site of the transfered hole 
ensures that the intermediate state Hamiltonian of the 
form Eq. (10 1 gives the correct scattering amplitude. To 
lowest order in tpd we hence find 



UdjUd + Up)^''{Up + 2A) 
2{Ud - ^)tld{2Ud + Up + 2A){Up + A) 



-4, (12) 



which results in r/ — —1.3 -again restricting ourselves 
to superexchange of the Anderson and Geertsma type. 
We see that to lowest order, the core-hole potential Uc 
does not appear in the well-screened intermediate state. 
From these microscopic considerations we conclude that 
the intermediate state Hamiltonian Eq. ( 10 1 is the correct 



one and higher order corrections to it are small because 
for the cuprates ?7 is a number of order unity. 

In a previous theoretical treatment we have shown in 
detail how to derive the cross section for RlXS-processes 
with a local core-hole using the UCL expansion^. For an 
incoming/outgoing photon with momentum qin/qout and 
energy ujf^/co^^^, we obtained t he cross se ction through 
the Kramers-Heisenberg relatiorPSEHSHMl as a function 



of energy loss ui 

q = qout - qin: 



and momentum transfer 



(Pa 



(Eldijj 



res \ 



Afi = Wrcs E 



Lo ~ LOf i) ) , with (13) 

/ rp 

{f\D\n) {n\b\t) 

Win - En - iT 



(14) 



The initial state \i) with energy Ei (which is used as 
reference energy: Ei — 0) is photo-excited to an interme- 
diate state which is described by the dipole operator D. 
The system can evolve through the intermediate states 
\n) with energy En (measured with respect to the reso- 
nance energy Wrcs) and, after the decay of the core-hole, 
end up in a final state |/) with energy Ef. Because the 
life time of the core-hole is ultrashort, we introduce an 
energy broadening F for the intermediate state. The de- 
tuning of the incoming photon energy from the -RT-edge 



is given by Win = ^fn ~ '^res- Finally, the delta function 
in Eq. ( 13 1 imposes energy conservation: the energy gain 
of the system uj fi — E f — Ei must be equal to the energy 



loss of the photon uj 



If F > En we can 



expand the amplitude Afi in a, powerseries. We assume 
that the energy of the incoming photon is tuned to the 
resonance (wjn = 0): 



oo 

^^-3^E(Z^(/I^(^.nt)' 



(15) 



1=1 



Note that we left out the I — term because it only 
contributes to elastic scattering. The leading order non- 
vanishing term in the sum is Z = 1, since the core-hole 
broadening is quite large compared to J. At the cop- 
per K-edge is 2r « 1.5 eV according to Refs. I19I20I and 
2F w 3 eV for the closely related ions Mn and Ge ac- 
cording to Refs. I21I22I which in either case is large com- 
pared to J. As in the three-band model rj = —1.3/ — 0.3 
eV for the well-/poorly-screened intermediate state, the 
largest value we find is rjJ/T « —0.22. Note that the 
UCL expansion therefore converges very well -even faster 
for the poorly-screened state than for the well-screened 
state (where \t]\ is larger). It is possible to directly in- 
clude a number of terms with Z > 2 in the cross section 
by using the expansion 



E 

1=1 



(gjnt)' 
F' 



E 

1=1 



Hi 
F' 



(16) 

with H' = vE^,JS,slJ,JS^■SJ. Since [Ho,D] - and 
^^0 I*) =0, all terms with Hq on the right can be safely 
neglected. Using Eq. (16), Afi simphfies to 
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with the scattering operator 



E s^"^ Jij Si • Sj 



(17) 



(18) 



From this equation we can deduce two important fea- 
tures. Firstly, indirect RIXS probes a momentum de- 
pendent four-spin correlation functior^. Secondly, Oq 
commutes with the z-component of total spin Sz , so the 
allowed scattering processes should leave Sz unchanged. 
Only an even number of magnons can be created or an- 
nihilated. 

To bosonize Eq. (|l8|, we split Oq in four parts: 



o.= E' 

i,j<£A 



-E' 



- E ■ 



- E •■■ (19) 

ieB, jeA 



Next, we rewrite this expression using LSWT as intro- 
duced in section In] Fourier transforming the result gives 
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' jAB 
, k+q/2 
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Jo-- -J. 



AA 

q 



^k- 



q/2) (ftk+q. 



Jo" 



/2"-k+q/2 



'^k-q/2^k+q/2 



^k-q/2"k+q/2 , 



'^k-q/2^-k-q/2 



(20) 



and we can write Oq in terms of the magnon operators 
using the inverses of Eqs. ([2| and (|3|. This leads to 



Oq = 0« 



0(2) 



(21) 



■(12). 

where Oq is a lengthy expression that contains the 
one/two-magnon scattering part. The next section deals 

with the two-magnon part Oq where two magnons 

are created or annihilated. The one-magnon part Oq^'' 
(where the change in the number of magnons is zero) is 
treated in section HVl 



III. TWO-MAGNON SCATTERING AT T = K 



At T = K, the system is in its groundstate, where 
no magnons are present: |i) = |0). Adding conservation 
of Sz, the only allowed scattering processes are the ones 
in which two magnons are created, so we consider the 



two-magnon part of the scattering operator of Eq. (21 1 
with 5=1/2: 



keA/BZ 



-J, 



jAA 



k+q/2 -^1^41/2 ^ Jq"^ ~ + + '^q j ("k+q/2l'k-q/2 + Uk-q/2Wk+q/2) + 

('^k+q/2 + -^k-q/i) (wk+q/2Mk-q/2 + l'k+q/2l'k-q/2) (ak+q/2/5-k+q/2 + "Lq/2/^lk-q/2) (^2) 



The two-magnon spectrum is shown in Fig.jsja). Several 
remarkable features can be seen. 

First of all the spectral weight vanishes at q ~ (0, 0) 
and q = (tTjTt), as can be seen in Fig. |4|b). This is in 
agreement with experimental observation^^^. The van- 
ishingof the RIXS intensity at q = is obvious: from 
Eq. ( 18 1 we see that at q = 0, Oq reduces to 2Hq (the 



factor of 2 arises from the fact that the sum in Eq. ( 18 ) 
is over all i and j). At zero temperature, \i) = |0) and 
consequently iJp |*) = -the RIXS intensity vanishes. 
At nonzero temperatures, Hq \i) = Ei\i) and according 
to Eq. (17) only elastic scattering occurs. It is easy to 



show that at q = (tt, tt) the RIXS intensity always van- 
ishes, regardless of the temperature or the form of Jij 
(as long as there is antiferromagnetic order). This holds 
because q = (tt, tt) is a reciprocal magnetic lattice vector: 



gjq-R. = 1 if R. is in sublattice A and e'^J'^* 



-1 if Hi 



is in sublattice B (assuming that at = (0, 0) we are 
in sublattice A). We obtain 

Oq— (tt.tt) — ^ ^ JjjSi'Sj — ^ ^ JijSj-Sj. (23) 

ieA,j leB.j 

Adding all terms where j G _B in the first term and j ^ A 
in the latter, we get zero. What remains is 

Oq=(7r,7r) = E Jij^i'^] ^ E '^ij^i'^j- (24) 



These terms cancel when applied to an initial state which 
is symmetric under the interchange of the sublattices. 

The other remarkable feature of the magnetic RIXS 
spectrum is its strong dispersion. This is apparent from 
Fig. |3ja) and|4]ja), showing the first moment (average 
peak position) of the spectrum. The calculations for 
the nearest neighbor Heisenberg antiferromagnet (see the 
dashed line in Fig.|4|a)) show that the magnetic scatter- 
ing disperses from about a; « around (0, 0) to w « 4 J at 
(tt, 0) and {tt/2, tt/2). Longer range couplings tend to re- 
duce (increase) the first moment of the RIXS spectrum if 
they weaken (reinforce) the antiferromagnetic order (see 
the solid line in Fig. 4|^a)). The observed dispersion in 
Fig. [sja) has a two- fold origin. It is in part due to the q- 
dependence of the two-magnon density of states (DOS), 
combined with the scattering matrix elements that tend 
to pronounce the low energy tails of the two-magnon 
DOS. In Fig. |3|b), it looks as if the two-magnon DOS 
has two branches. The most energetic one around q = 
is strongly suppressed by the matrix elements throughout 
the Brillouin zone (BZ). 

The consistency at q = (0, 0) and q — (tt, tt) of the 
theoretical results and experimental data was already no- 
ticed, but at other wave-vectors the agreement stands out 
even more. The data on La2Cu04 for q = (tt, 0) shows 
a peak at around 500 meV, precisely where we find it on 
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FIG. 3: RIXS spectrum (a) and two-magnon DOS (b) for a 
nearest neighbor Heisenberg antiferromagnet with exchange 
interaction J as a function of transferred momentum q for a 
cut through the Brillouin zone (c). The dashed line indicates 
the magnetic BZ boundary. 



the basis of a nearest neighbor Heisenberg model with 
J = 146 meV - a value found by the analysis of neu- 
tron scattering datsF^. Similar agreement is found at 
q = (O.GtTjO) and q — (0.67r, O.Gtt).^'' Even better agree- 
ment is found when we take into account the second and 
third nearest neighbors and ring exchange according to 
the neutron data. The ring exchange interaction, which 
we treat on a mean field level, simply renormalizes first- 
and second-nearest neighbors exchang^^I! 




(0,0) (Jt,0) (71,71) (0,0) 

(b) 

FIG. 4: First moment (a) and total spectral weight (b) of the 
RIXS spectrum. The solid lines are obtained by using interac- 
tion strengths determined from neutron data (next neighbor 
coupling J — 146.3 meV, second and third neighbor couplings 
J' = J" = 2 meV and ring exchange Jc = 61 meV).'^^ The 
dashed lines have only nearest neighbor interaction. 



In Fig. [5] we compare the results for the two-magnon 
scattering intensity with experimental data,G2! using the 
interaction strengths determined from neutron dataP^l, 
for three values of q in the BZ. Note that we use the 
wave- vector independent renormalization factor Zc here, 
that takes into account some of the magnon-magnon in- 
teractions.'^ This simply changes the energy scale by a 
factor Zc ~ 1.18 but does not affect the intensity of the 
spectrum. Each panel shows the theoretical prediction 
(dashed line), the theory convoluted with the current in- 
strumental resolution (solid line), and the experimental 
data. The only free parameter in the theoretical spectra 
is the over-all scale of the scattering intensity. We find it 
to vary by a factor of 2.5 comparing different q's, which 
is within the error bars of the experimeniPs!. 

Many qualitative features such as the occurrence of in- 
tense peaks at the magnetic BZ boundary and the large 
dispersion characterizing the total spectrum are in ac- 
cordance with our earlier result and the results of Na- 
gao and IgarashP^. The spectra of Ref. UHl taking two- 
magnon interactions partially into account, show slight 
quantitative differences with respect to our results: the 
RIXS peaks soften and broaden somewhat as a conse- 
quence of the magnon-magnon interaction, particularly 
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FIG. 5: RIXS intensity for various points in the BZ. Each 
figure contains the bare theoretical data (dashed line), the 
convolution with experimental resolution (solid line), and 
the experimental data from Ref. [I7l For these figures, we 
used J = f46.3 meV, second and third neighbor couplings 
J' = J" — 2 meV and ring exchange Jc = 61 meV. The latter 
contribution is evaluated theoretically using a mean field ap- 
proximation. These values were found in neutron scattering 
experiments;^^ These experiments were analyzed using the 
wave- vector independent renormalization factor Zc = 1.18, 
which is also used to generate the theoretical curves. The 
theoretical intensity is scaled independently in each figure to 
match the experiment. The overall scale factors differ at most 
by a factor 2.5, which is cornparable to experimental uncer- 
tainty in absolute intensities.'^Sl 



for the (tTjO) point. The range of the dispersion in the 
spectrum is therefore smaller (the mean uj/J varies be- 
tween 1 and 3 instead of 1 and 4) . 



IV. FINITE T: SINGLE-MAGNON 
SCATTERING 



The S^gf symmetry allows scattering processes where 
no additional magnons are created. In the finite tem- 
perature case, an initial magnon of momentum k can be 
scattered to k— q. The one-magnon part of the scattering 
operator, within LSWT, takes the following form: 



keMBZ 



jAB ^ jAB _ jAA _ jAA ^ jAA ^ jAA^~^ (^^^^_^ ^ 



tAB 



(25) 



We choose to concentrate on the basic case where the 
only non-vanishing interaction is the nearest-neighbors 
coupling J, for a 2D Heisenberg antiferromagnet with 
5"= 1/2. 

In the low temperature regime, a few magnons of low 
momentum k are present in the system. Their energy 
can be approximated for T — > by letting k — > 0: ~ 
•\/2J|k|. In this limit Mk and Wk can be substituted by 
the following approximate expressions: 



In order to calculate the one-magnon contribution to 
the cross section, we have to evaluate the scattering 
amplitude expressed by Eq. (17 1. In the low temper- 



ature case we can consider a one-magnon initial state 
|i) — a]^ |0).'^The only contribution to A^^^ comes from 
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the final state with a single magnon of momentum k — q 



-{Jk + ./k-q)(wkWk-q 
g 

I j= (Jo + Jq) ("q - 

V2V2 



- WkWk-q) 



(27) 



where we used the condition |k| ^ |q| and inserted the 
expressions of Eqs. ( 26 1 for Uk and Wk, retaining the lead- 



ing order term in |k|. 

These approximations allow the analytic evaluation of 
the scattering intensity. At finite T, the cross section is 
given by 



dClduj 



1 



.(1) 



5{uj~Ef + Ei). (28) 



For k w 0, and by taking the continuum limit, we obtain 



dflduj 



cx F(q) / dk 
Jmbz 



.dk 



y ^fie^ _ ^ '^(^-ek-q + ek), 
(29) 

where we defined P(q) = S'^ ( Jq + Jq) (itq — I'q)^- In 
the low temperature limit, the Bose factor goes to zero 
rapidly for high |k|, so the only substantial contribution 
to the integral comes from |k| w 0. Therefore we can 
extend the domain of integration to the entire fc-space. 
Replacing ek with its approximate expression in the limit 
of low |k|, and assuming polar coordinates, we obtain 



cx P(q) / dk- 
Jo t 



fc2 



1 



/2Jfc) (30) 



Eq, which 



Note that we used the replacement ek-q 
breaks down at q = and the BZ corners. This integral 
can simply evaluated to be 



d^duj 

and the spectral weight for T/ J <C 1 is 



(eq - ^) 



— du oc P(q) 



(31) 



(32) 



The T"^ behavior also shows up in the numerical eval- 
uation of Wi (without assuming |k| ^ |q|), as shown 
in Fig. [6] as a function of the transferred momentum q, 
for various temperatures (dashed lines). According to 
the considerations discussed in the previous section, the 
RIXS intensity is vanishing for (7r,7r). The average peak 
position and the peak width are expected to be modified 
as a function of temperature. We can easily estimate 
these modifications by evaluating the first moment 



30C(3) 



(33) 
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FIG. 6: Comparison between spectral weight for single- 
rnagnon scattering Wi (dashed lines) for various temperatures 
and zero temperature two-magnon scattering W2 (solid line), 
all obtained numerically. The figure displays the behavior 
from Eq. (32 1 for the single-magnon intensity. For La2Cu04 
J ~ 146 meV, and at room temperature we have Pit J ~ 5.8. 



and the variance 



) (X 1 . 



(34) 



We conclude that the peak position is shifted from eq 
towards lower values of lo, by an amount that grows lin- 
early with T and at the same time the peak broadens 
proportional to T. 

We now determine the relative intensity of the one- and 
two-magnon scattering processes. Even if a direct com- 
parison is not possible, since the one-magnon and the 
two-magnon peaks occur at different lost energies w, it is 
useful to compare the one-magnon and the two-magnon 
total spectral weight for the 2D Heisenberg antiferromag- 
net. The latter is evaluated numerically at T — 0, and 
the former at various temperatures without making the 
approximation k — q w — q. In Fig. [6] we plot the two- 
magnon (solid line) and the one-magnon weight for differ- 
ent temperatures (dashed lines). At room temperature, 
the one-magnon weight is one or two orders of magni- 
tude smaller for almost every value of q and is expected 
to decrease with decreasing T, according to Eq. ( 32 1 . 



This allows us to conclude that the two-magnon scatter- 
ing is the dominant process at low temperatures. A rough 
estimate for the temperature at which the one-magnon 
process becomes significant gives a value of ~ 1 eV in the 
case of La2Cu04, which is well above room temperature. 
These results support the conclusion that two-magnon 
scattering dominates the magnetic RIXS intensities in 
this material observed by J. P. Hill and coworkers^^. In 
other materials this of course needs not necessarily be so, 
depending on the temperature at which the experiments 
are performed. One can expect for instance interesting 
RIXS scattering signals from high temperature param- 
agnons. 




LEADING CORRECTION TO ULTRASHORT 
LIFETIME APPROXIMATION 



The ultrashort core-hole lifetime (UCL) expansion 
offers a systematic way of calculating the Kramers- 
Heisenberg relation Eq. (13 1. In this section we calcu- 



late the leading correction term to the two-magnon cross 
section in the UCL approximation. This is especially rel- 
evant at q = (0, 0) where the intensity is vanishing to first 
order, but non-zero to second order. The leading order 
correction is taken into account by including all terms 
up to ©((TyJ/r)^) in Eq. (16). Again we can include a 



number of extra correction terms by using an expansion 
of the type 



1=1 



1=1 



E 

1=2 



(35) 



The contribution of the last term to the UCL scattering 
amplitude is 



iT + LO 



(36) 

This scattering amplitude that corresponds to this term 
is non-zero at q = 0, which can be easily checked in linear 
spin-wave theory. The reason is that the resulting scat- 
tering operator at zero transferred momentum does not 
commute with the Hamiltonian. For the LSW analysis 
we make use of the identity 



- jE4S«-S, +const. (37) 



We drop the constant because it does not contribute to 
inelastic scattering. For simplicity, we only take near- 
est neighbor interactions into account. The last term in 
Eq. ( 37 1 is proportional to the first order result for the 



scattering amplitude, which has already been analyzed 
in LSWT. The other term can be treated in LSWT too, 
and yields a two-magnon contribution to the scattering 
amplitude at zero temperature of: 



4r2 iT - 



E(/i/(k.q)^ 



("k^k+q + Mk+qi'k)"],;/?! 



k-q I 



(38) 



with /(k,q) — —6{cosqx+cosqy)+4:COskxCOs{ky + qy) + 
4:C0skyC0s{kx + Qx) + 2cos(2A:a; + q^) + 2cos{2ky + qy). 
Since the phase of the first order amplitude differs from 
the second order amplitude by 7r/2, there is no interfer- 
ence of these terms. The consequence is that the lead- 
ing corrections to the first order scattering intensity are 
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FIG. 7: The leading order correction to the scattering amph- 
tude does not interfere with the first order. Fig, (a) shows 
the contribution to the cross section from Eq. ( |36| |. The full, 
corrected cross section is shown in Fig. (b). There is an ap- 
preciable correction only at q = 0. 



down by a factor {rjJ/T)^ « 0.06 for the well-screened 
intermediate state. This makes the ultrashort core-hole 
lifetime approximation a viable way of computing mag- 
netic RIXS spectra. The contribution Eq. ( 36 1 is shown 
in Fig.[7]^a), and the full cross section in Fig. [7|jb). Only 
at q = there is an appreciable difference from the first 
order result shown in Fig. [3] (a). At q = (tt, tt), there is 
again no intensity, which can be understood by the same 
argument as for the first order result in section [llTl 



VI. CONCLUSIONS 

We derived the two-magnon scattering cross section 
which is measured in magnetic RIXS at the Cu if-edge, 
taking advantage of a series expansion in the ultrashort 
core-hole lifetime (UCL) of the intermediate state. In 
the context of LSWT, we calculated the magnetic RIXS 
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spectrum for a 2D S' = 1/2 Heisenberg antiferromagnet, 
in the more general case where the superexchange is not 
Umited to nearest neighbors. Our results strongly suggest 
a multi-magnon scattering scenario, where two-magnon 
excitations are created in the system as a consequence 
of the modifications in the superexchange interaction in- 
duced by the core-hole potential. 

Our results for the two-magnon scattering agree very 
well with experimental data on La2Cu04. The vanishing 
of the RIXS intensity for the elastic case q = (0, 0) and 
the antiferromagnetic point q = (tt, tt) is recovered. The 
latter feature turns out to be a consequence of an under- 
lying symmetry property of the scattering operator and 
does not depend on the range of the exchange interac- 
tion. The excellent quantitative agreement between our 
results and experiments is testified by the occurrence of 
an intense peak at q = (tt, 0) for uj sa 500 meV. We have 
generalized the theory to include also finite-temperature 
scattering, for which we find that also one-magnon pro- 
cesses contribute. For La2Cu04 at room temperature the 
single magnon spectral weight is very small compared to 
two-magnon scattering. 

The subleading order in the UCL expansion of the 
cross section is shown to be of order ©((t] J/F)^) smaller 



than the first order result. This makes the UCL approx- 
imation a rigorous method for this case to calculate the 
Kramers-Heisenberg relation. The introduction of longer 
range interactions (according to data from neutron ex- 
periments) improves the correspondence between theory 
and magnetic RIXS experiments on La2Cu04. The gen- 
eralization of the analysis to doped systems will be an 
interesting next step towards understanding multi-spin 
correlations in the spin liquid phase of the high-Tc su- 
perconductors. 
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